Dynamic canonical suppression of strangeness in transport models 
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It is investigated whether canonical suppression associated with the exact conservation of an 
U(l)-charge can be reproduced correctly by current transport models. Therefore a pion-gas having 
a volume-limited production and annihilation cross section for nn «-» KK is simulated within two 
different transport prescriptions for realizing the inelastic collisions. It is found that both models 
can indeed dynamically account for the canonical suppression in the yields of rare strange particles. 
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I. INTRODUCTION 

As the properties of hadronic gases provide important 
theoretical means of describing heavy ion collisions, there 
exist several models to address this problem. Among 
the most prominent are statistical thermal models for 
equilibrated hadronic gases. 

One result is a suppression in the yields for rare par- 
ticles when treating the conservation of a corresponding 
U(l)-charge exactly [l|,|2>|3|. The importance of such a 
treatment is fortified by recent experimental results 
The suppression effect can eventually be formulated in 
terms of a volume-dependence . 

A dynamical interpretation has recently been of- 
fered by the formulation and solution of kinetic master- 
equations, which are capable of covering the time evolu- 
tion of the system as well |2| . 

In the following we will demonstrate that standard 
hadronic transport models dynamically show canonical 
suppression. More explicitly the U(l)-charge considered 
is the strangeness whose net value is taken to be zero 
throughout our calculations. Due to simplicity we con- 
sider only inelastic reactions of the type irir <-> KK, 
where the kaons and anti-kaons bear strangeness +1 and 
— 1 respectively. 

An abundant number of kaons justifies a grand canon- 
ical description of the system, in which the strangeness 
is conserved only on the average (Nk) — (N K ) = 0. The 
density of particles having strangeness +1, the kaons, can 
then be computed 
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where Z K denotes the relativistic one-particle partion 
function for non-interacting particles (kaons in this case) 
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A rare number of kaons however demands the conserva- 
tion of strangeness to be treated exactly, i.e. Nk — Nk = 
0. A consistent realization of this constraint eventually 
leads to a canonical partition function describing the sys- 
tem 0, Q and thus to the following kaon density 
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with the canonical suppression factor < ry < 1 being 
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and x — 2Z K . Pursuing the kinetic master equation 
approach yields the same results for the equilibrated kaon 
densities 0. 

Investigating the canonical suppression factor Q more 
closely, using the asymptotic behaviours 
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one finds that in the grand canonical limit the kaon den- 
sity n,K is independent of the reaction volume, whereas in 
the canonical regime, with the number of kaons (Nk) <C 
1, it scales linearly with the volume as x oc V. This 
volume dependent behaviour of the kaon density pro- 
vides a convenient reference for comparison between sim- 
ulation and theory. 



II. DYNAMICAL SIMULATIONS 

The simulation setup consists of a large box of 20 fm 
side length holding a relativistic gas of pions. This pion- 
gas provides a heat bath for a much smaller reaction vol- 
ume of variable size centered within the large box. In- 
side this likewise box-shaped reaction volume processes 
7T7T «-> KK are allowed, covering all possible isospin states 
of pions and kaons. The kaons are reflected by the walls 
of the reaction volume, whereas these walls are perme- 
able for the pions. After equilibration, the kaon density 
within the reaction volume should be governed by 
which holds as a reference for the transport model results. 
Therefore different sizes of the inner reaction volume are 
simulated and the number of kaons (Nk) is extracted by 
averaging over many timesteps. The timesteps are suffi- 
ciently large such that no correlations are to be expected. 

We implemented this scenario using two different types 
of transport descriptions - the microscopic transport 
model UrQMD Q and a realization of a stochastic trans- 
port model borrowed from a recent and elaborated parton 
cascade 6] . The former model makes use of a geometrical 
interpretation of cross sections in order to solve transport 
equations, whereas the latter relies on the calculation of 
transition probabilities. 
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The UrQMD model provides full space time dynam- 
ics for hadrons and strings. It is a non-equilibrium 
model based on the covariant propagation of hadrons 
and strings and the geometric interpretation of cross sec- 
tions. For our studies it was modified such that only 
reactions tttt <-> KK together with all possible elastic 
2 <-> 2 processes remained possible. We had to choose 
sufficiently small cross sections such that the mean free 
path remained large compared to the interaction length 
A m .f.p. = (71022) _1 3> v c r 22/ 7r - Otherwise the difference 
in the collision times of the involved particles viewed from 
the lab frame causes a decrease in the collision rates as 
pointed out in |(|. As different densities are involved, for- 
ward and reverse reactions are not affected in the same 
way and the change in reaction rates leads to a shifted 
equilibrium value of the kaon density. For details see • 

The simulation of 2 <-> 2 processes within the stochas- 
tic method is based on the calculation of a collision prob- 
ability for each possible particle pair per unit volume Ax 3 
and unit time At via 

At , , 

P22 =v re ia 2 2-^. (6) 

v re i denotes the relative velocity and 022 is the cross sec- 
tion for the considered 2^2 process. Any so obtained 
probability is then compared with a random number be- 
tween and 1 to decide whether the collision should take 
place or not. 

The initial conditions in any model are chosen such 
that the pion gas acquires a temperature of T = 
170 McV. The appropriate number of pions and the to- 
tal energy of the system are calculated via the use of a 
grand canonical partition function for pions alone. Our 
heat bath volume of 8000 fm corresponds to a popula- 
tion of 1348 pions bearing a total energy of 747.5 GeV. 
Initially each pion is assigned the same fraction of the 
total energy, giving one half of the particles momenta 
in the positive x-direction, while the remaining particles 
start out bearing momenta in the negative cc-direction. 
The spatial distribution is random. 

A fundamental and necessary verification of the simu- 
lations' reliability is of course a check for thcrmalization. 
Figure^riicely demonstrates that the energy distribution 
of the pions becomes thermal in both models. 

Figure |2 then displays our results in terms of the 
extracted kaon density plotted versus the reaction vol- 
ume along with the theoretical predictions for the grand 
canonical (JIJ and canonical behaviour. The minor 
fluctuations in the results provide a visual indication for 
the irrelevance of statistical errors. It is manifest that 
canonical suppression is reproduced by both models, i.e., 
as expected the kaon yield is suppressed with respect to 
the grand canonical limit for small reaction volumes and 
thus small numbers of produced kaons. 

The canonical suppression for the number of kaons 
(Nk) being considerably smaller than 1, as pointed out 
in [2] , originates dynamically from an enhancement of the 
annihilation process by 1/ (Nk) as compared to standard, 
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FIG. 1: Logarithmic energy spectra for the simulated pion- 
gases. 
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FIG. 2: Kaons density versus reaction volume as extracted 
from simulations with UrQMD (diamonds) and the stochas- 
tic cascade (triangles). For comparison the dashed line indi- 
cates the grand canonical behaviour. The solid line shows the 
canonical volume dependence of the kaon density as expected 
from © for a temperature of T = 170 MeV. 
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grand canonical formulation of the Boltzmann equation. 
The reason is that any kaon in the system requires the 
existence of a corresponding anti-kaon due to strangeness 
conservation. Thus the probability of finding a particle 
anti-particle pair turns into a highly correlated condi- 



tional probability. The so enhanced annihilation proba- 
bility then leads to the canonical suppression in the kaon 
yields. As shown by our results, this behaviour is auto- 
matically included in the considered transport models. 
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